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ABSTRACT 



We investigate which polynomials can possibly occur as fac- 
tors in the denominators of rational solutions of a given par- 
tial linear difference equation (PLDE). Two kinds of polyno- 
mials are to be distinguished, we call them periodic and ape- 
riodic. The main result is a generalization of a well-known 
denominator bounding technique for univariate equations to 
PLDEs. This generalization is able to find all the aperiodic 
factors of the denominators for a given PLDE. 

Categories and Subject Descriptors 

1.1.2 [Computing Methodologies]: Symbolic and Alge- 
braic Manipulation — Algorithms 

General Terms 

Algorithms 

Keywords 

Difference Equations, Rational Solutions 

1. INTRODUCTION 

Several algorithms in symbolic computation depend on a 
subroutine for finding the rational solutions of an ordinary 
linear difference (or differential) equation (OLDE), and sev- 
eral algorithms are known for implementing such a subrou- 
tine [1, 2, 4, 11, 13, 14, 6, 8, 9]. On a conceptual level, the 
typical approach for finding rational solutions can be divided 
into three steps. In the first step, one constructs a polyno- 
mial Q such that the denominator q of any potential solution 
p/q must divide Q. This polynomial Q is called universal de- 
nominator or denominator bound. In the second step, the 



*Supported by the Austrian Science Fund (FWF) grants 
P20347-N18 and Y464-N18. 

'Supported by the Austrian Science Fund (FWF) grants 
P20162-N18 and P20347-N18. 



Permission to make digital or hard copies of all or part of this work for 
personal or classroom use is granted without fee provided that copies are 
not made or distributed for profit or commercial advantage and that copies 
bear this notice and the full citation on the first page. To copy otherwise, to 
republish, to post on servers or to redistribute to lists, requires prior specific 
perinission and/or a fee. 

ISSAC2010, 25-28 July 2010, Munich, Germany. 
Copyright 2010 ACM 978-1-4503-0150-3/10/0007 ...$5.00. 



Carsten Schneider 

RISC 

Johannes Kepler University 

4040 Linz (Austria) 

cschneid@risc.jku. at 

universal denominator is used to transform the given equa- 
tion into a new equation such that P is a polynomial solution 
of the new equation if and only if P/Q is a rational solution 
of the original one. In the third and final step, the polyno- 
mial solutions P of the transformed equation are determined. 

The first algorithm for computing universal denominators 
in the case of OLDEs with polynonrial coefficients was pro- 
posed in 1971 by Abramov [1] (see Section 2 below for a sum- 
mary). It has been generalized to g-difference equations [3], 
to matrix equations [5], and also to equations whose coef- 
ficients belong to domains other than polynomials. For ex- 
ample, Bronstein [7] and Schneider [12] have observed that 
a universal denominator can be constructed also when the 
coefficient domains are difference fields which can be used 
for representing nested sums and products (IIE-fields). For 
such domains, the situation is more involved. There is a 
need to distinguish between "normal" factors of the univer- 
sal denominator which can be found very much like in the 
usual polynomial case, and "special" factors which have to 
be constructed by some other means. 

In the present article, we consider partial (i.e., multivari- 
ate) linear difference equations with polynomial coefficients 
(PLDEs). Our ultimate goal is the construction of a univer- 
sal denominator for potential rational solutions of a given 
PLDE. Like in the univariate case with sophisticated coef- 
ficient domains, there are two kinds of factors to be distin- 
guished. As a matter of fact, some parts of the denominator 
cannot be bound at all. For example, the equation 

/(?H-l,fc) = /(n,fc+l) 

has {n-\-k)~°' as a rational solution, for any a G N, and there 
is obviously no finite polynomial Q that would be a multiple 
of (n -I- k)" for all a G N. We will call factors that may 
exhibit such "special" behaviour periodic. Our main result 
is that we can construct for any given PLDE a polynomial 
d such that every aperiodic factor of any potential solution 
p/q must divide d. 

Such a bound on the aperiodic factors of the denominators 
does not directly give rise to a full algorithm for finding ra- 
tional solutions of PLDEs, but it can be considered as a step 
in this direction. For a full algorithm, besides of the bound- 
ing of the periodic parts of the denominator, also the entire 
question of how to find (all) polynomial solutions of a PLDE 
in the third step is wide open and far from being settled. But 
even if these parts have to remain open for now, our aperi- 
odic denominator bound is useful in practice. When it comes 
to solving an actual equation, possible periodic factors in a 



solution can often be guessed by inspection, their multiplici- 
ties can be determined by trial and error, and degree bounds 
for polynomial solutions can be established heuristically. A 
reasonably tight universal denominator, on the other hand, 
cannot be as easily obtained on heuristic grounds. 

2. THE UNIVARIATE CASE 

Before entering the multivariate setting, let us summa- 
rize Abramov's classical denominator bound for univariate 
equations. We will introduce on the fly some notions and 
notations needed later. 

Let K be a field of characteristic zero and let K[n] and 
K(n) denote the ring of univariate polynomials and the field 
of rational functions in n with coefficients in K, respectively. 
Write N for the shift operator acting on K[?i] and K(n) via 

Nq{n) :=g(n+l). 

The objects of interest are difference equations of the form 

aoy + aiNy-\ + amN'^y = f (1) 

where ao, . . . , am, f £ K[n] (ao, a-m / 0) are given and y £ 
K(n) is unknown. 

The denominator bounding problem is as follows: given 
ao, . . . ,am, f G KM, find Q G K[n] \ {0} such that the 
denominator of any solution y G K(n) of (1) divides Q. 
Abramov's denominator bounding algorithm [1] is an effi- 
cient way of computing 



a divisor of Ili^o N''am, and therefore the denominator of y 
must be a divisor of 



gcd(n^'fflo,n^"'""'«"0, 



where 



s :— max 



{i>0:gcd(iV'ao,iV"'"a„)/l} 



(2) 



(3) 



is the dispersion of ao and N~"^a,n- It is efficient in the sense 
that the gcd is constructed without explicitly calculating the 
products. 

To see that this bound is correct, write (1) in the form 



N"^y=±(f-J2<^^N\_ 



Shifting this equation by s gives 



N'^+'y ■ 



N'>a„ 



^(^V-E(^'"')(^'+^y) 



By repeatedly using the recurrence, the terms N^^"y appear- 
ing on the right hand side can be reduced to smaller shifts 
of y so that for certain polynomials b,bo, . . . , b,n~i we have 

At this point we rely on the following result. 

Theorem 1. [1] For any solution y = ^ £ K(n) of (1), 

max{i > : gcA{q, N"q) / 1} < s. 

This theorem ensures that the denominator of a solution 
y cannot contain two factors u,v with u — N"^'''v, and this 
in turn implies that no denominator of any of the N^y on 
the right can have a common factor with the denominator 
of N'^^'y. Therefore the denominator of iV'+^y niust be 



n^v^ 



niv 



(4) 



By an analogous argument, now rewriting y in terms of 
higher shifts, it can be shown that the denominator of y 
must divide 11^=0 N^ag. As both bounds must hold simulta- 
neously, the correctness of Abramov's bound is established. 
The argument given here is not exactly Abramov's origi- 
nal one, but follows a proof which to our knowledge first 
appeared in [5]. The equivalence of the two approaches is 
shown in [8]. 

3. THE MULTIVARIATE CASE 

For the rest of this article we will be concerned with adapt- 
ing the univariate reasoning sketched in the previous sec- 
tion to the multivariate setting. From now on, we consider 
polynomials and rational functions in r variables ni , . . . , n^. 
Wherever it seems appropriate, we will use multiindex nota- 
tion, writing for instance n for ni, . . . ,nr, etc. 

We define shift operators Ni, . . . , Nr acting on K[n] and 
K(n) in the obvious way: 

iVig(ni,.. .,nr) 

= q{ni,. . . ,ni-^i,ni + l,ni+i,. . . ,nr). 

For i — {ii, . . . ,ir) £ Z^ we will abbreviate 

N'q~Nl'N^^---K^q. 

A partial linear difference equation (PLDE) is an equation 
of the form 



s6S 



(5) 



where S C Z'" is finite and nonempty (called the support or 
the shift set or the structure set of the equation), / G K[n] 
and fls G K[n] \ {0} (s G S) are explicitly given polynomi- 
als, and y G IK(n) is an unknown rational function. The 
polynomial as is called the coefficient of A'''^ (or simply of s). 
In the following section, we generalize the notion of dis- 
persion to multivariate polynomials, and we will define the 
notions of periodic and aperiodic polynomials. After that, 
in Section 5, we will show how to predict all the aperiodic 
factors of a rational solution y of (5). 

4. SPREAD AND DISPERSION 

The notions of spread and dispersion are related to the 
question whether two polynomials can be mapped to one 
another by a shift. In the multivariate case, we now allow 
independent shifts in all directions. 

Given two polynomials p,q G K[n], we say that they are 
shift equivalent if there exists i = (ii , . . . , v) G Z'' such that 
p(n) — q{n + i). In operator notation, p and q are shift 
equivalent iS p = N'q. 

Definition 1. Let p,q £K[n]. The set 
Spread(p,(7) := { (ii, . . . , v) G Z'' : 

gcd(p,iV^.--iV;-g)/l} 
is called the spread of p and q. The number 

Dispj.(p, q) := max{ \ik\ : (Ji, ■ ■ ■ , ir) G Spread(p, q) } 



IS called the dispersion of p and q w.r.t. k £ {1, . . . , r}, and 

Disp(p, q) := max(Dispi {p,q),..., Disp^(p, q)) 

is called dispersion of p and q. (By convention, max A := 
— oo if A is empty and max^ := cxj if A is unbounded.) 

The polynomial p is called periodic if Spread(p, p) is infi- 
nite and aperiodic otherwise. 

Note that this definition does not exactly correspond to 
the definition stated before for the univariate case. While 
there, the definition depends on whether one shifts to the left 
or to the right [1, 7], our definition takes all directions into 
account. This makes the reasoning below a little simpler. 

Example 1. Let 

p={n + k") {{n + if + (k- 3f) {k-n + 3), 
g = ((n + 2f + {k- if) {{n - 2f + {k + 7f) {2k - 3n) . 
Then, by inspection, 

Spread(p,Q) = {(2,-l), (-2,7), (1,2), (-3,10)}, 
Disp(p, q) = 10. 

Both p and q are aperiodic. An example for a periodic poly- 
nomial is n — k, because, again by inspection, 

Spread(?i — fc, ji — fc) — {{i,i) : i £ Z } 
= {...,(-l,-l),(0,0),(l,l),(2,2),...} 

is infinite. 

In the univariate case, the spread of two polynomials can 
be found as the set of all integer roots of the polynomial 
res„(p(n),g(n + z)) G K[i]. This is no longer possible in the 
multivariate setting, for in the case of several variables, com- 
mon roots no longer correspond to common factors. Never- 
theless, it turns out that the multivariate spread as defined 
as above can be effectively computed. Let us consider the 
somewhat simpler situation of irreducible polynomials first. 
In this situation the spread cannot take on any cardinality: 

Lemma 1. Let p,q £ IK[n] be irreducible and aperiodic. 
Then | Spread(p, g)| < 1. 

Proof. Suppose p and q are such that | Spread(p, q)\ > 1. 
Then there exist two different multiindices (ii, . . . , v) and 
(ji, •■•,>) with 

gcd(p, Nl' • • • iV^g) / 1 and gcd(p, iVf • ■ • AT^g) / 1. 

As p and q are irreducible and irreducibility is preserved 
under the shifts n^ i— >■ n^ + 1, we have in fact 

cp = iVJi ■■■N^Sq^ iVf • • • iV^'-g 

for some c G K \ {0}, hence 

q = Nl^-'^ ■■■N'/~'''q, 

in contradiction to the assumption that g is aperiodic. D 

It is not essential for the lemma that we consider only 
shifts of integer distance. More generally, if p and g are two 
irreducible polynomials, then, by the same argument, the 
number of vectors i £ K'^ with p(n) = g(n -I- i) is either 
or 1 or infinite (K refers to the algebraic closure of K). 



These vectors i £ K can be found by making a brute force 
ansatz. For a variable vector i = (ii, . . . , v), force 

p(n)-g(n + i)=0 
and compare coefficients with respect to n to obtain an alge- 



braic system of equations in ii. 



over K. The solutions 



will form an affine linear space over K, for whenever i, j £ K'^ 
are such that p(n) — g(n + i) and p(n) = g(n + j), then 

p(n) — g(n + i -I- a(i — j)) for all a £ Z, 

and since the solution set is Zariski-closed, what is true for 
all a £ Z must also be true for all a £ K. 

The spread of two irreducible polynomials can therefore 
be computed by first determining a basis of the affine lin- 
ear space of all possible shifts i £ K*^ mapping one given 
polynomial to the other. By taking the radical to remove 
nontrivial multiplicities, it is ensured that a basis of the lin- 
ear space can be read off from a Grobner basis. In a second 
step, we filter out from this affine space the vectors which 
have integral coordinates only: 

Algorithm 1. /npui; p, g £ K[n] irreducible 
Output: Spread(p, g) 

1 S := Coefficients(p(n) - g(n -h i), {n}) C K[i] 

2 G := Gr6bnerBasis(Radical(S),degrevlex(i)) 

3 Choose a basis B of the (Q vector space generated by 

the coefficients of the elements of G, say B = 
{6i,...,64CK. 

4 For each g £ G, let g^^\ . . . ,p'''' £ (Q[i] be such that 

g = 6ig'^' -f &25*^^ H h &dg''''. (Note: At this 

point all g £ G are linear, and so are all the g^ ' .) 

5 S:=[j{g^'\g^'\... ^^''^^ 

geG 

6 return kera(5') 



'} 



Note that the algorithm avoids the need of solving systems 
of diophantine equations by exploiting a priori knowledge on 
the structure of the solution set. 

The case of non-irreducible polynomials is easily reduced 
to the former case by considering all pairs of factors. To be 
precise, let p, g £ K[n] \ {0} be any polynomials, and let 



(ll 12 ■■■Is 



P = Pl P2 ■■■Pr , q 

be their factorization into irreducible factors. Then 

r s 

Spread(p,g) = [J IJ Spread (p», g^). 

1=1 j=i 

In short, given p,q £ K[n] \{0} we can compute Spread(p, g) 
and therefore also Disp^(p, g) and Disp(p, g). 

Every given polynomial p can be split uniquely (up to 
constant multiples) into a factorization p = uv where u is 
periodic and v is aperiodic. We call u and v the periodic 
and aperiodic part of p, respectively. As we can factor poly- 
nomials and compute, as described above, their spread, this 
decomposition can be computed. 

Example 2. 1. p = n + k, q — n + 2k. Here we have 
p{n, k) - q{n -\-i,k+ j) = (-i + 2j) - k 
and G — {1}, hence Spread(p, g) — 0. 



2. p = n+ V^k, q = n + y^k + 3 - 2y^. Here we have 

p{n, k) - q{n + i,k+ j) = -i - V^j - 3 + 2\/2 

and G = {-i-^/2j- 3 + 2^/2}. WithB = {1,V^} we 
get S = {— i — 3, —j + 2}, from which we obtain 

Spread(p,Q) = {(-3,2)}. 

3. p = 3fc^ + 6fcn - 7fc + 3n^ - 7n + 1, q = 3/c^ + 6fc7i - 
13fc + 3n^ — 13n + 11. Here we have 

p{n, k) -q{n + i,k + j) ^ -{i + J - l)(3i + 3j - 10) 
- 6(i + j- l)k + 6{i +j- l)n 

and G = {i+j - 1}. With B = {1}, we get S = G 
from which we obtain 



Spread(p,g)= (q) +'z(\ 



5. APERIODIC FACTORS IN DENOMINA- 
TORS OF SOLUTIONS 

In this section we solve the following problem. Given 
a nonempty finite shift set S C Z' with coefficients as £ 
K[n] \ {0} (s G 5"), find a polynomial d G K[n] \ {0} such 
that for any solution y = -^ of (5) with p,q,u £ K[n] where 
q is aperiodic and u is periodic, we have 

q\d. 

Such a d is called an aperiodic universal denominator (or 
aperiodic denominator bound) of (5). 

First we generalize Theorem 1, i.e., for any solution y G 
K(n) of (5) we bound the dispersion of the aperiodic de- 
nominator part of y. To be more precise, we first bound the 
dispersion w.r.t. one component Ui in n. 

Lemma 2. Let S '^ HT be finite and nonempty and let 
Us G K[n] \ {0} for s G 5", and f G K[n]; let a'^ be the 
aperiodic part of a^. Let i G {1, . . . , r}. Define 

k := max{|ai - bi\ : (ai, . . . ,ar), (6i . . . ,br) G S} 

and let 

A —{(si, . . . , Sr) : (si, . . . ,Sr) e S and 

3(tl,. . . ,tr) e S S.t. ti- Si = k}, 

B ={(si, . . . , Sr) : (si, . . . , Sr) G 5 and 
3(ii, . . . ,tj.) £ S s.t. Si — ti = k}. 

Define 

Si := max{Dispj(as, Afj at) : s G A and t G B}. 

Then for any solution y ~ -^ & K(n) of (5) with periodic 
part u and aperiodic part q we have 

Dispjg) =Disp.(q,g) < Si. 



Proof. Suppose that d :— Disp;(g) > Si. Then we find 
irreducible factors u and v of q such that 






N'^u — V 



(6) 



for some e = (ei, . . . , e,) G Z'' with d = d. Consider all the 
factors N^u and N^v occurring in q where the ith entries in 
u and V are 0. We choose now those factors from q where u 
and V are maximal w.r.t. lexicographic order; these factors 
are denoted by u' and v' , respectively. 

First suppose that u' divides one of the polynomials as with 
s G ^. As 5 is not empty, B is not empty. Therefore, we 
can choose that polynomial aw with w — {wi, . . . ,Wr) G 
B such that {wi, . . . ,Wi-i,Wi+i, . . . ,Wr) is maximal w.r.t. 
lexicographic order. By (5) we can write 



iV' 



y^ 



aw V 



ses\{w} 



(7) 



Now observe that the factor N^v' does not occur in the 
denominator of any N^y with s G S* \ {w}: if N^v' oc- 
curred in such a denominator, s G -B by construction (recall 
that u and v' have maximal distance d in the i-coordinate 
among all factors in q and that v' with N'^v' is shifted 
maximally by k among all possible choices from S in di- 
rection of the i-coordinate since w G -B; so only if s G -B 
is necessary to guarantee that N^v' is a factor of N^q and 
hence of the denominator of N^y). But then by the as- 
sumption that (uii, . . . , Wi-i, Wi+i, . . . , Wr) is maximal w.r.t. 
lexicographic order (among all possible choices) and that 
v' = N^v with V = {vi, . . . ,Vi-i,0,Vi+i, . . . ,Vr) is maximal 
w.r.t. lexicographical order, it follows that N^v' can only 
occur in the denominator of N^y. Summarizing, the factor 
N^v' does not occur in the denominators of N'^y for any 
s £ S \ {w}, and since /, Os G K[n], the common denom- 
inator of / — X]ses\|w> o.s^^y does not contain the factor 
N^v' . Moreover, since u is a factor of as for some s G A, 
and since w G -B and d > Si, N^v' cannot be a factor of 
as for any s G -B. In particular, our aw is not divisible by 
N^v' . Overall, the common denominator of the right hand 
side of (7) cannot contain the factor N^v' which implies that 
the denominator of N^y is not divisible by N^v' . Thus the 
denominator of y, in particular q is not divisible by t;'; a 
contradiction. 

Conversely, suppose that u' does not divide any of the 
polynomials as with s G A. Then by similar arguments as 
above (the roles of A and B exchanged), we derive again a 
contradiction. Therefore Si < Dispj((jr). D 

Example 3. Ln the generic univariate case (1) (r = i = 
1) the shift set is S = {0, 1, . . . , m} C Z^ and for the sets A 
and B from Lemma 2 we have A = {0} and B — {m}. In 
this particular instance, Lemma 2 corresponds to Theorem 1. 

A bound of the dispersion for the multivariate case is given 
in the next theorem. 

Theorem 2. Let S '^HT be finite and nonempty and let 
as G K[n] \ {0} for s G S, and f G K[n]. Then one can 
compute an s G N U {— cxj} with the following property: For 
any solution y — -^ £ K(n) of (5) with periodic part u and 
aperiodic part q we have 

Disp(g) = Disp(g, q) < s. 



(8) 



Proof. Compute the values Si for i £ {1, . . . ,r} as de- 
scribed in Lemma 2; the spread can be computed with Al- 
gorithm 1. By taking s — max(si, . . . , s^) the property 
Disp(q) = Disp(g, q) < s is guaranteed. D 



In order to derive an aperiodic denominator bound, we 
adapt the idea presented in Section 2. Namely, we will 
choose an appropriate point p £ S and express A''''^ for 
any solution y £ K(n) of (5) in terms of N^y for points 
s £ S' which are sufficiently far away from p. To be more 
precise, for any s > we can explicitly write 



N^y: 






(9) 



for some polynomials b, bi £ K[n] and for finite sets W, S' C 
"E^ with the following property: the distance of the points S' 
to p is at least s. Then by taking s as in Theorem 2, we will 
be able to conclude that Higw-p -^''^p is an aperiodic de- 
nominator bound of N^y, and consequently riicw N'~^''ap 
is an aperiodic denominator bound of y. 

Such appropriate points p from S can be chosen as follows. 
Let 5 C Z"^ be a finite set. A point p £ S is called a 
corner point (or extreme point) of S, if there exists an affine 
hyperplane H (codimension 1) which contains p and where 
all other points S \ {p} are situated in one of the two open 
half spaces determined by the hyperplane. 




Such an affine hyperplane H is called border plane of S for 
p, and a vector being orthogonal to H and directing to the 
half space of the points S \ {p} is called inner vector. 

Note that the corner points are the extreme points of the 
convex hull generated by S, and they can be computed by 
simple linear algebra; for further details see, e.g., [10]. 



Example 4. In the generic univariate case (1) (r — 1 
and S = {0, 1, . . . , m}) the corner points are and m, and 
the border planes are {0} and {m} with inner vectors 1 and 
— 1, respectively. 

More generally, if we are given a finite set 5 C Z*" with 
(0, . . . , 0) £ S and max{di : (di, . . . , dr) £ S'} > for each 
1 <i <r, there are at least r -\-l corner points. 



For our denominator bound construction we start with the 
following simple lemma. 



Lemma 3. Let S (^ Z^ be a nonempty finite set and let 
p £ S be a comer point together with a border plane H and 
inner vector v. Consider any hyperplane H' which is parallel 
toH. Then for any p' £ H' [MT i/ie poinis S'-|-(p'-p)\{p'} 
are all outside of H' in the half space determined by the 
direction of v. 




Proof. Since H + (j>' — p) = H', H' is a border plane of 
S + (p' — p) for p'. This proves the lemma. D 

By iterative application of the previous lemma we obtain 
the following theorem. 

Theorem 3. Let S ^71/' be nonempty and finite, and let 
Os £ K[n] \ {0} for s £ S, and f £ K[n]. Let p be a comer 
point of S with a border plane H and inner vector v. Then 
for every s > there exist finite sets 



H^ C Z" n [J {H + ew) and 

0<e<s 

S'cZ'-n |J(// + (s + e)v), 



(10) 

(11) 



and polynomials b,bi £ K[n] such that for any solution y £ 
K(n) of (5) the relation (9) holds. 




Proof. We show the theorem for a generic solution y £ 
K(n) of (5). Hence the ingredients b, fe, £ K[n] and W, S' un- 
der consideration will hold for any specific solution as stated 
in the theorem. For S :— S \ {p} we have 



N^y: 



/-Es6S«=^=2/ 



by (5). If 5* = {}, take S" = {} and W = {p} and we are 
done. Otherwise, let p' £ S be such that the distance to H is 
minimal. DeRne H' := H +{p' -p). li H' <Z {H + {s + e)v : 
e > 0}, we are again done with W = {p} and S' — S \ {p}. 
Now suppose that H' C {H + ev : < e < s}, and let 
{pi , . . . , Pfc } = H' n S (by construction p' £ H' n S) and 
define Si = 5* \ {pi, . . . , p^}. Then we can write 



N^y: 
By (5) 



/ - Eses, «siV=j/ - E.e{p„....p,> «=iV=y 



ATP 



APl-POn 



(12) 



(13) 



In particular, by Lemma 3 the points from S-f (pi — p) are all 
outside of H' in the half space determined by the direction 



of V. Thus we substitute (13) into (12) and can express N^y 
in the form 



/' - E.es; <N^y - E 



SG{P2, ■■-,?/,.} 



a.N^'y 



N^y: 

for some /' £ K[n] and a'^ e K[n] with S[ = 5iU(S+pi-p) 
After k — 1 further reductions, we end up at the form 



N^y = 






for some 6,6s £ K[n] and with Wi = {p. Pi, • • ■ ,Pfc} and 
S'l = S'iUUi<i<fc 'S'+(pi— p). Note that all points pi,. . . , Pfe 
which are closest to H have been eliminated. Now we re- 
peat the construction from above until we enter in the base 
case given in the beginning of the proof. This completes the 
proof, n 

Note that all the ingredients W, S' , b £ K[n] and the 6s for 
s £ W can be computed explicitly. However, for getting an 
aperiodic denominator bound, we only need W. The proof 
of Theorem .3 delivers the following simple algorithm. 

Algorithm 2. Input: A finite nonempty set S C Z'' and 
a corner point p of S with a border plane H and inner vector 
v; s > 0. 

Output; A finite set W^ C Z' with (10) such that there are 
S' with (11) and 6, 6i £ lK[n] such that (9) holds for any 
solution y G K(n) of (5). 

1 S-5\{p};S'-S\{p};W^:={p} 

2 while IJ (-^ + ""^^ n SV {} do 

0<e<s 

3 Let {pi, . . . , pk} be the points m S' which have min- 

imal distance to H . 

4 W^:=WU{pi,...,Pk} 

5 5':=S'\{pi,...,Pk} 

6 S*' := S' U U 5" + (pi - p) 

l<i<fe 

7 enddo 

8 return W 

Finally, we end up at the following theorem which tells us 
how we can compute an aperiodic denominator bound. 

Theorem 4. Let S '^HT be finite and nonempty, and let 
as e K[n] \ {0} for s e S, and f G K[n]; let a'^ be the 
aperiodic part of a^. 

Take s G N U {— oo} s.t. for any solution y ~ -^ & K(n) 
of (5) with periodic part u and aperiodic part q we have (8) . 
Let p be a comer point of S with a border plane H and inner 
vector V with |v| > 1. Let W be the output of Algorithm 2 
with input H , p and s. Then 



n 



iV^a' 



sGW-2p 

is an aperiodic universal denominator of (5). 

Proof. Let y = -^ & K(n) be a solution of (5) with pe- 
riodic part u and aperiodic part q. By construction of W it 
follows that there are 5*' and 6, 6i G K[n] with (11) and (9). 
Since |v| > 1, the distance of the points S' to p is larger than 
s. Observe that N^q must occur in the denominator of the 
right hand side of (9). Using (8), the aperiodic denominator 



parts in N^q with s G 5*' cannot contribute to the aperi- 
odic denominator part of N^q. Hence only the polynomial 
risew- ^"'^p is responsible for N^q, i.e.. 



Thus ri: 



sew-2p 



sG W — p 

N^Op is a universal denominator of (5). D 



Example 5. In the generic univariate case (1) (r = 1 
and S = {0, . . . , m}) the coefficients are given by Oi for < 
i < m. First, we take the corner point m with the border 
plane H = {m} and inner vector —1. Note that with (3) it 
follows that (8) holds for any solution y = - £ lK(n); see 
also Example 3. Applying Algorithm 2 with the input S, m, 
H , —1, we get W = {tti, m — 1, . . . , m — s}. Hence Theorem 4 
delivers the universal denominator bound (here we have only 
aperiodic factors) 



n 



N'a„ 



s6M''-2m 

which agrees with (4). Similarly, if we take the corner point 
with the inner vector 1, we get W = {0, 1, . . . , s} and obtain 
the universal denominator bound 

s 

n N'ao = l[N*ao. 

seW-2'O i=0 

Combining these two estimates produces (2). 

From the point of view of application the following re- 
marks are in place. 

1. s G NU {— oo} can be computed by Theorem 2 and by 
applying Algorithm 1. 

2. Applying Algorithm 2 we can compute the finite set 
W (^ W,^; here we remark that different choices of the 
border plane H might lead to sets W of different size. 
Exploiting the particular structure of S gives room for 
improvement. 

3. Suppose that we are given k corner points with corre- 
sponding border planes. Then by Theorem 4 we end 
up at different aperiodic universal denominators, say 
di, . . . ,dk G K[n]. Then taking 

gcd(di,...,dfc) 

leads to a sharper universal bound. 

4. The coefficients as with s G S are often available in 
factorized form. Then also the d^s are obtained in fac- 
torized form, and the gcd-computations boil down to 
comparisons of these factors and bookkeeping of their 
multiplicities. 

Combining the aperiodic denominator bounds for different 
corner points gives the following result. 

Theorem 5. Let S (^ Z^ be finite and nonempty and let 
as G K[n] \ {0} for s G S, and f G K[n]. iei pi, . . . ,pfc G 
S be corner points of S. If the denominator of a ratio- 
nal solution of (5) contains an aperiodic irreducible factor, 
then shift equivalent factors occur m each of the coefficients 



Proof. By Theorem 4 aperiodic denominator bounds can 
be derived by the corner points pj in the form Yli -^'op, , 
respectively. Hence an aperiodic denominator bound of (5) 
can be written in the form 

d=gcd(niv'ap,,...,n^'«pj- 



If a rational solution contains an aperiodic irreducible factor 
h, then h is also contained in d. Hence /i or a shift equivalent 
factor occurs in each of the a„. , . . . , a„,. . D 



The following special cases are immediate. 

Corollary 1. Let S C Z' be fimte and nonempty, and 
let as G K[n] \ {0} for s € S, and f G K[n]. Let s,t € S 
be two comer points and let a^ and a^ be the aperiodic parts 
of the coefficients as and at, respectively. If'Disp{a's,a[) — 
— oo, then the aperiodic denominator part of any rational 
solution of (5) is 1. 

Corollary 2. Let S C Z'' be finite and nonempty, and 
let as G K[n] \ {0} for s G S, and f G K[n]. // there is a 
corner point of S whose coefficient has no aperiodic factor, 
then the aperiodic denominator part of any rational solution 
of (5) IS 1. 

Besides these structural consequences, Theorem 5 pro- 
vides the following improvement of our aperiodic denomi- 
nator bound algorithm. To be more precise, Lemma 2 and 
thus Theorem 2 can be improved in the following way. In 
the proof of Lemma 2 we assume that there are irreducible 
factors u and v in the denominator of the solution y G K(n) 
of (5) such that (6) for some e = (ei,...er) G "Z^ with 
Ci = d where d is larger than Si. By the choice of Si this 
leads to a contradiction. Now we exploit in addition Theo- 
rem 5: the factors it and v can be only factors that occur -up 
to shift equivalence- in each coefficient of the corner points 
pi , . . . , Pfc . Hence it suffices to choose Si as summarized in 
the following proposition. 

Proposition 1. Let S C S'' be finite and nonempty, and 
let Os G K[n] \ {0} for s G 5, and f G K[n]. Let pi, . . . ,pfe 
be corner points of S, and let a^ be the aperiodic part of as 
whose factors are present -up to shift equivalence- in each 
coefficient of the comer points. Let i G {1, . . . , r}. Define 

k := max{|ai — &;] : (ai, . . . , Or), (&i . . . , br) G S} 

and let 

A ={(si, . . . , s,.) : (si, . . . , Sr) G 5* and 
3{ti,...,tr) G 5" s.t. ti - Si = k} 

B ={(si, . . . , s,.) : (si, . . . , Sr) G 5* and 
3{ti,...,tr) G 5* s.t. Si -ti =k}. 

Define 

Si ;= max{Disp^(as, A''" a^) : s G A and t G B}. 

Then for any solution y = -2- G lK(n) of (5) with periodic 
part u and aperiodic part q we have 

Disp,(g) =Disp.(q,g) < Si. 



6. EXAMPLES 

Example 6. Consider the recurrence 
{2kn + l)(6fc^ + 12fc - W - 4?i + 5)/(n, k) 

+ {2kn + 4k + l)(6fc^ + lOfc + A^n -h 8n - 7)/(n -f 1, k) 

- (2kn + Sn + l)(6fc^ + 24k + 4n - 20n - 7)/(n, A: -f 2) 

- {2kn + 4fc -I- 8n -f 17) 

X (6fc^ + 22fc - 4n 4- 16n + 45)/(n -f 1, fc -f 2) = 0. 

The maximum spread among the coefficients of this recur- 
rence is s = 4. 

Every point in the shift set {(0, 0), (1, 0), (0, 2), (1, 2)} qual- 
ifies as a corner point. We choose p = (0, 0) as corner point 
and let H be the plane through p orthogonal to v = (1, 1). 

Algorithm 2 delivers the set 

{(0, 0), (1, 0), (2, 0), (3, 0), (4, 0), (5, 0), (6, 0), (7, 0), (8, 0), 
(0,2), (1,2), (2, 2), (3, 2), (4, 2), (5, 2), (6, 2) 
(0,4), (1,4), (2, 4), (3, 4), (4, 4), 
(0,6), (1,6), (2, 6), 
(0,8)} 
OS W, from which by Theorem 4 it follows that 

8 4-ri/21 

Y[ n N'K^^2kn+l){6k^ + 12k-'in^ -in + b) 

i = j=0 

is a universal aperiodic denominator. 

Taking instead (1,2) as corner point gives the aperiodic 
denominator bound 



Yl n N'-^ K^^~'^{{2kn-\- 4k + 8n -It 17) 

X (6A:^ -f 22fc - 4n -\- 16n + 45)) . 
The greatest common divisor of the two polynomials is 
{2kn + l)(2(fc + 2)n + l)(2fc(n + 1) + l)(2(fc + 2)(n + 1) + 1). 
This is exactly the denominator of the actual solution 

3fc-|-n 

(2fcn+l)(2(fc + 2)n + l)(2fc(n + l) + l)(2(fc + 2)(n + l) + l) 

of the recurrence. 

The computation could have been simplified by disregarding 
the factors p = 6fe^ -|- 12fc - 4n^ - 4n + 5 and g = 6fc^ -I- 22fc - 
4n^ -|- 16n -I- 45. Because of 

Spread(p,g) = {}, 

they cannot contribute to the universal denominator (com- 
pare Theorem 5). 

Example 7. Some comer points may be easier to handle 
than others. As an example, consider the equation 

{k^ + n + l){2k^ + 4fc^ + Ak^n^ + Ak^n + 6fc^ + Akn^ 

+ 8kn + 9fc -f 2n* -f 4n^ + 6n^ -f 3n -I- 4)f{n, k) 

- {k^ + 4k + n + 5)(2fc^ -I- 4k^ + Ak^n + Ak^n + Qk^ - 1 
+ Akn -2kn-k + 2n + An^ + &n ~ 2n)f{n, k + 1) 

- (2fc + l){n + l)(fc^ +n^ -hAn + 5)f{n + 1, fc) = 0. 

Without any computation, it can be deduced that any poten- 
tial aperiodic factor in a denominator must be a shifted copy 



o/ fc^ + 71^ + 4n + 5 = (71 + 2)^ + A:^ + 1. Indeed, a rational 
solution of the equation is given by 

fc^+n^ 

(fc2 + n2 + 1) {{k + 1)2 + n2 + 1) (fc2 + (n + 1)2 + 1) ■ 

Example 8. Theorem 4 ^^ not sufficient for predicting pe- 
riodic factors of a denominator. As an example, consider the 
equation 

2{k + n + l)/(n, k)-{k + 3n + 8)/(n, fc + 1) 
- (5fc + 3n + 12) f{n + l,k)+3ik + n + 5) fin + 1, fc + 1) 
+ (fc + n + 5)/(n + 2,fc) =0. 

Possible choices for corner points are (0,0), (0, 1), (1, 1) and 
(2,0). Because of 

Spread(fc + n + 1, 3n + fc + 8) = {}, 

one might be tempted to believe that only trivial denomina- 
tors can occur in a solution. However, the equation admits 
the nontrivial rational function solution 

n^ +k'^ 



{k + n + l){k + n + 2)(fc + n + 3) ' 

Observe that although not every corner point contains a 
shifted copy of k -\- n -\- 1, there is still some comer point 
which does. 

This need not be the case, as indicated by the example 

/(n+l,fc)-/(n,fc + l)=0 

already mentioned in the introduction. This example, how- 
ever, is special because the shift set {(0, 1), (1, 0)} belongs 
to a proper affine subspace of TIT , and whenever this is the 
case, say the shift set belongs to a subspace L C Z^ , then 
for every vector v = [vi, . . . ,Vr) € L^ \ {0} the polynomial 
p = wiWi + • ■ • + VrUr has the property that f is a rational 
solution of the equation if and only if fp" is, for any a £%. 
In particular, if there exists a rational solution at all, then 
there also exists one whose denominator does not contain p. 
Apart from this exceptional situation, we observed on all 
the examples we considered that periodic factors of a denom- 
inator appeared (possibly as a shifted copy) in at least one of 
the coefficients corresponding to some comer point of the 
shift set. This at least suggests the periodic factors of these 
coefficients as plausible guesses for the periodic part of the 
denominator bound. 

7. CONCLUSION 

There are polynomials in several variables which may have 
an infinite spread. Such polynomials are called periodic, 
while polynomials that cannot have infinite spread are called 
aperiodic. Each polynomial can be split into a periodic and 
an aperiodic part. 

We have shown that for partial linear difference equations 
with polynomial coefficients it is possible to determine all 
the factors that may possibly occur in the aperiodic part of 
the denominator of a rational function solution. The con- 
struction is a generalization of the corresponding result for 
univariate equations. It probably admits further generaliza- 
tion to (jf-equations or equations whose coefficients belong to 
a IlE-field. As it stands, Theorem 4 will tend to produce 
only a rough bound for the aperiodic part of the denomina- 
tor, but we have pointed out several refinements for improv- 
ing the efficiency of the computation on concrete examples. 



Eventually, it would be interesting to see whether it is possi- 
ble to come up with an Abramov-style algorithm for directly 
computing the greatest common divisor of all the individual 
bounds obtained from each corner point. Until now, we have 
not succeeded in constructing such an algorithm. 

It also remains open how to bound periodic factors of the 
denominator. The situation illustrated in Example 8, which 
seems to be typical both for random equations and for equa- 
tions coming from applications, indicates that our result for 
aperiodic factors does not directly extend to periodic ones. 
On the other hand, it also indicates that an equation typi- 
cally provides some hints for the periodic factors in the de- 
nominators of its rational solutions. This is useful for making 
plausible heuristic guesses. It also gives some hope that at 
least for certain types of equations the periodic part of a de- 
nominator can be found algorithmically. This needs further 
investigation. 
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